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I. INTRODUCTION 



Studying the structure of the nucleon and its excita- 
tions is essential in the development of an understanding 
of the strong interaction at the confinement scale [Ij . On 
a fundamental level there exists only a limited under- 
standing of the relationship between Quantum- Chromo- 
Dynamics (QCD), the field theory of the strong interac- 
tion, and models that use effective degrees of freedom, 
such as the large variety of constituent quark models or 
alternative hadron models [3, [j, 0, la, Q ■ Experimentally, 
we still lack sufficiently complete data sets that can be 
used in phcnomcnological analyses to unravel the internal 
structure of the nucleon's excited states by determining 
their electromagnetic transition amplitudes for different 
distance scales. Precise cross section and polarization 
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data for ir° production off protons, which arc needed 
to study the transition from the nucleon ground state 
to the A(1232), have become available in recent years. 
For a recent review see Ref. J7|- These data have been 
used to extract the electromagnetic transition multipoles 
for this state over a wide range of photon virtuality Q 2 
@, @, E3, EH d 111 EI El 111- These accurate data 
have been confronted with recent results from quenched 
Lattice QCD calculations [l7l ll8J| . 

The A(1232)P33 excitation corresponds to the transi- 
tion between isospin 1/2 and isospin 3/2 ground states 
of the non-strange baryons and is characterized by a spin 
flip in the ground state. To study the radial and or- 
bital excitations of the nucleon, we need to measure the 
transition form factors for higher mass states. In the 
mass region above the A(1232)P33, there is a cluster 
of resonances around 1.5 GeV that consists of 3 isospin 
1/2 nucleon excitations. The 7V(1440)Pn, often called 
the Roper resonance, is thought of as the first radial 
excitation of the nucleon. The two other states, the 
A r (1520).Di3 and AT(1535)5n, arc orbital excitations of 
the nucleon. By measuring the transition form factors of 
these states, we probe the radial and orbital wave func- 
tion of the nucleon that describe fundamental properties 
of the nucleon structure at the constituent quark and 
meson-nuclcon scales. In the mass range around 1.7 GeV, 
there is a cluster of nine resonances that consists of a 
number of isospin 1/2 and isospin 3/2 states in various 
configurations of radial and orbital excitations of the 3- 
quark system. These resonances overlap with each other 
and contain higher partial waves. This mass range is 
also characterized by significant non-resonant contribu- 
tions to the production cross section. In order to sep- 



arate these states, and isolate them from non-resonant 
contributions, more experimental information is required 
than in the case of the A(1232). Measurements of dif- 
ferent isospin channels, e.g. pir° and nir + final states, 
are essential. Measurements with polarized beam and 
polarized target, and measurements of the proton recoil 
polarization, arc highly sensitive to interferences of res- 
onant and non-resonant amplitudes contributing to pion 
production. These data arc indispensable for isolating 
individual resonances and determining their excitation 
strengths, within the framework of any model. 

The mass range above the A(1232) has recently been 
systematically mapped out in measurements of unpo- 
larized cross sections of the electroproduction processes 



ep -> enn+ [19|, |2{fl , ep -> epn |9|, ep -> eprj [21|, [22J, \2c 
and ep — ► epTT + ir~ [2J, |25[, and measurements of beam 
spin asymmetries using polarized electron beams [19 . |26J ] . 
Detailed measurements with polarized beams and polar- 
ized targets were carried out in inclusive electron scat- 
tering in the resonance region j27|, [28|, [29| , and a limited 
number of data points were measured for the ep — > enn + 
reaction with polarized proton targets [30(. However, 
no data are available for the processes ep — > epir and 
ep — > epir . The intent of the current research has been 
to fill this gap and to provide precise and detailed data 
from measurements for the pir° final state using both po- 
larized electrons and a polarized proton target, and to 
systematically map the mass region from W = 1.1 GeV 
to W = 1.7 GeV and in the range 0.187 < Q 2 < 0.770 
GeV 2 . 




FIG. 1: (Color online) Schematic diagram of 7r°-nucleon 
electro-production, e represents the incident polarized elec- 
tron, el is the outgoing electron, 7* is the virtual photon, 
and p and p are the nucleon in the initial and final states, 
respectively. 



where d£t e = sin9 e d9 e d(j) e is the electron solid angle and 
dCl* = sin 9*d6*d<f>* is the solid angle of the meson in the 
center of mass. Here 



II. FORMALISM 

The ep — ► epir reaction is schematically shown in Fig- 
ure [TJ The incident polarized electron is given by the 4- 
vector p e = (p e ,Ei), the outgoing electron is emitted with 
angles <p ei 6 e and 4- vector p' e = (p e ',Ef) 7 and the nu- 
cleon initial and final states are given by p v = (0, M) and 
Pp = (Pp'Ep), respectively. The virtual photon 7* with 
4- vector q = (q,u>) where q = p e — p e ' and u> = Ei — Ef , 
is characterized by the squared four momentum transfer: 



Q 2 = -q 2 = -{uo 2 -\q | 2 ) = 4EjEi sin 2 -± (1) 
and the degree of transverse polarization: 



e = ll+2 ^! tail ^- 



(2) 



The invariant mass of the hadronic system is given by: 



W = be +Pp -Pel = V M2 + 2MW ■ 



(3) 



In terms of these variables the differential cross section 
can be written as: 
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denotes the 'photon equivalent energy' necessary for a 
real photon to excite a hadronic system with center-of- 
mass (cm.) energy W, and a is the fine structure con- 
stant. 

The differential cross section for pion production by a 
virtual photon da/dfl* can be written as a sum of four 
terms: 
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where k is the momentum of the pion, h is the electron 
hclicity, and P is the target proton polarization. The 
first term doo/dn* represents the unpolarized cross sec- 
tion, while the remaining terms dcr e /dn*, dot/dn* , and 
dcr e t/dn* arise when beam, target, or both beam and 
target are polarized, respectively. Here 
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is the 'real photon equivalent energy' in the cm. frame. 
These cross sections can be written in terms of response 



J 



functions R using the formalism of reference [3l| as: 
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(10) 



is the frame-dependent longitudinal polarization of the 
virtual photon. The d 1 is the angle between the direc- 
tions of the target polarization and the virtual photon. 
The asymmetries are then defined as: 



A , --1± 

-tt-et — • 

00 



(11) 



where o§ = dao/dQ*, a e = da e /d£l* , a t = da t /dQ*^ and 
a et = da et /dQ*. 

If integrated over the azimuthal center-of-mass angle 
(f>* , the double spin asymmetry A et can be parameterized 
as [H: 



VT^ 
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where r\ = tan0 7 y / 2e/(l + e) and R is the longitudinal- 
transverse cross section ratio ox/ ' &t The structure func- 
tion A\ is the virtual photon helicity asymmetry, 
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while A 2 is a longitudinal-transverse interference term. 

III. REACTION MODELS 

In recent years, several models have been developed 
that aim at accurately reproducing the experimental 



data. In section IVII we compare some of our data with 
calculations based upon model descriptions such as the 
Dubna-Mainz-Taipei(DMT) model 0], a recent version 
of the MAID model 0, and the Sato-Lee model (SL) 0]. 
In addition, a unitary isobar model (UIM) was developed 
by the Yerevan- JLab group [5|, l|| that contains many 
features of MAID, but incorporates different energy- 
dependences of the background amplitudes at higher en- 
ergies. This approach allows the fitting of experimental 
cross sections and polarization asymmetries and the ex- 
traction of resonance transition amplitudes from a large 
body of data. We briefly summarize the main features 
of these models. They are discussed in more detail in 
Ref. 0. 

MAID and related models are based on an isobar de- 
scription of the single pion production process. They 
incorporate non-resonant amplitudes described by tree- 
level Born terms, and also include p and to i-channel pro- 
cesses that are relevant mostly in the region of higher 
resonances. Figure [5] shows the diagrams contribut- 
ing to the reaction ep — > epir at low and intermedi- 
ate energies. The vertex functions for the virtual pho- 
ton coupling to hadrons are parameterized according to 
their respective on-shell form factors for which there is 
prior experimental information. Resonances are parame- 
terized by a phenomenological description using a rel- 
ativistic Breit-Wigner form with an energy-dependent 
width. The total amplitude for single pion production 
is unitarized in a K-matrix formulation. Only single 
channels are included, and multi-channel effects such as 
77V — ► (pN, 7rA) — + ttN, which could be important in the 
second and third resonance regions, are neglected. From 
an experimental viewpoint, the attractive feature of these 
descriptions is flexibility that allows adjusting parame- 
ters such as electromagnetic transition form factors and 
hadronic couplings as new experimental information be- 
comes available. However, all of these descriptions lack 
significant predictive power, and a comparison with new 



data will tell us more about how well electromagnetic 
and hadronic couplings have been parameterized, rather 
than about the intrinsic structure of the nuclcon. 



(a) 



(b) 





(c)-£ 




FIG. 2: Tree-level diagrams of single piv electroproduction 
for (a) s— channel resonance production, (b) s— channel nu- 
cleon exchange and, (c) t— channel u meson exchange. 



Dynamical models, such as the SL and DMT models 
start from a consistent Hamiltonian formulation. In these 
models, the non-resonant interaction modifies the reso- 
nant amplitude. The SL model provides the most con- 
sistent description of the interaction but is currently lim- 
ited to the region of the A(1232) resonance, while in the 
DMT model the resonance amplitudes are parameterized 
according to a specific Brcit-Wigner form that simplifies 
the inclusion of higher resonances. The s-channel reso- 
nance parameterization in the DMT model is similar to 
what is used in the isobar descriptions such as the MAID 
and UIM approaches. Once the transition form factors 
have been extracted from the data, their interpretation in 
terms of the intrinsic structure of the nucleon must then 
involve comparisons with nucleon structure models, such 
as the many versions of the Constituent Quark Model's 
(CQM), and Lattice QCD calculations. 



IV. EXPERIMENTAL SETUP 

Measurements were performed in 2000-2001 using the 
CEBAF Large Acceptance Spectrometer (CLAS) |33( at 
the Thomas Jefferson National Accelerator Facility (Jef- 
ferson Lab), located in Newport News, Virginia. The ex- 
periment consisted of a longitudinally polarized electron 
beam with energies from 1.6 GeV to 5.6 GeV and currents 
ranging from 1 to 5 nA, scattering off longitudinally po- 



larized protons and deuterons in ammonia targets, 15 NBi3 
and 15 ND3. The beam was rastered in a spiral pattern 
of 1-1.2 cm diameter over the target to avoid destroying 
the target polarization due to beam heating. The beam 
had an average polarization of 70%, and the polarization 
direction was reversed with a frequency of 30 Hz. 

The experiment used two frozen ammonia targets, 
15 NH 3 and 15 ND 3 , and 12 C, 4 He, and 15 N targets for 
background studies. The targets were held in a 1 
K liquid-helium cooling bath and a 5 T longitudinal 
magnetic field generated by a pair of superconducting 
Helmholtz coils, which blocked particles scattered at 
polar angles between 45° and 70°. The protons and 
deuterons in the ammonia were polarized via Dynamic 
Nuclear Polarization and the polarization achieved was 
about 75% for protons and 25% deuterons. 

The scattered particles were detected using the CLAS 
detector, described in great detail in Ref. [331. The CLAS 
detector consists of 6 superconducting coils that produce 
a toroidal field, around the beam direction, which bends 
charged particles in the polar direction but leaves them 
nearly unaffected in the azimuthal direction. The polar- 
ity of the magnetic field was set to bend negative particles 
toward the beam line. Each of 6 sectors, delimited by the 
coils, are equipped with three drift chamber (DC) regions 
to determine the trajectories of the charged particles, a 
Cerenkov counter (CC) for electron identification, scintil- 
lators counters (SC) to measure the time-of- flight (TOF), 
and an electromagnetic calorimeter (EC) used to identify 
electrons and neutral particles. The detector covers 80% 
of the azimuthal angle. The DC and SC subtend polar 
angles from 8° to 142°, while the CC and EC cover up 
to 45°. 



V. DATA REDUCTION AND ANALYSIS 

In this paper we report the measurement of target and 
double spin asymmetries for the p(e, e'p)iro channel using 
the 1.6 GeV beam on 15 NH3 target data set. 



A. Channel identification 

The exclusive final state p(e, e'p)ir° was selected by 
identifying an electron, a proton, and a missing pion. 
Electrons are identified at the trigger level by requiring 
a coincidence signal from the CC and the EC, and the 
offline track reconstruction reduces the electron candi- 
dates by geometrically matching signals from the DC, CC 
and EC. Electron identification was further improved in 
the data analysis in order to reduce contamination from 
negative particles such as 7r~ and other sources of back- 
ground. The z-vertex position was used to remove events 
originating from scattering off the target windows. In or- 
der to improve the resolution of the vertex position, the 
data were corrected to account for the effect on the ver- 
tex reconstruction due to the beam rastcring. Figure [3] 




CQ.1.2 



-60 -40 -20 

z-vertex [cm] 



FIG. 3: Number of electron events as a function of the vertex 
z-position. The dashed and the solid lines represent the dis- 
tribution before and after the raster correction, respectively. 
The vertical lines indicate the cut applied to the data to re- 
move the peaks coming from the scattering off the exit win- 
dows of the target. 



shows the vertex position before and after the correction, 
and the applied cut, —58 cm < z-vertex < — 52 cm, to 
select events originating inside the 15 NH3 target. 

To separate electrons from negative pions, a combina- 
tions of cuts on the Cerenkov and the EC signals were 
imposed. The energy released by electrons traversing the 
EC is proportional to their momentum, and about 30% 
of this energy deposition is measured by the scintillator 
sheets in the EC (sampling fraction). To reject negative 
pions, a cut on the total EC visible energy as a function 
of the momentum was applied. The energy deposited in 
the EC is also measured separately for the inner part 
and the outer part. A cut to reject events with a low en- 
ergy deposition in the inner part of the calorimeter was 
applied. This allows further identification of electrons, 
which shower mostly in the inner part of the EC, while 
pions are minimum ionizing and lose most of their energy 
in the outer part. 

In CLAS, positive hadrons are identified using the mo- 
mentum and the path length information from the DC 
and the timing information from the SC. Figure U shows 
the cut applied to select protons, which appear well sep- 
arated from positive pions at the energy of this experi- 
ment. The missing pions were selected by a 2ct cut on the 
missing mass M\ = \p e +p p — p e > —p p >\ 2 spectrum around 
the mass of the pion. In addition, to reduce multi-pion 
background, events with detected charged particles other 
than the electron and the proton in the final state were 
rejected. 
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FIG. 4: j3 vs p for positively charged particles. The lines show 
how protons are selected in the analysis. 



B. Fiducial cuts and acceptance corrections 

The objective of this analysis is to extract asymme- 
tries, for which the precise knowledge of the acceptance 
is not as crucial as for differential cross section measure- 
ments. Nevertheless, the acceptance can still affect the 
results. To calculate asymmetries, the ratio of counts 
integrated over ranges of cf>*, 9*, Q 2 , and W is formed, 
and since the acceptance depends upon these variables, it 
does not necessarily cancel. Fiducial cuts were applied to 
the data in order to remove regions where the efficiency 
varies by more than one order of magnitude, and the 
acceptance was calculated analytically for the geometri- 
cal regions denned by the fiducial cuts. This procedure 
fully corrects the data for the acceptance except for an 
overall efficiency factor that is assumed to be constant 
within the fiducial cuts, and which will therefore cancel 
out when taking the ratios in the asymmetries. 



C. Experimental definition of the asymmetries 

For each hclicity configuration, the number of events 
can be written as: 
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where the arrows in AT-j-p, N-^, N^, and N±± indicate the 
beam and the target polarizations, respectively. Since the 
target was the same for all configurations, the parame- 
ters D\ through D4 differ only by the total accumulated 
charge Q. P® and P 4 a are the beam and target polariza- 
tions for runs with a positively polarized target, and P b 



and P t b for negative target polarization runs. The beam 
polarizations for the groups (N^,N^) and (N^i,Nn) are 
assumed to be equal because the polarization of the beam 
was reversed every half second, and therefore each pair 
was taken at the same time. On the other hand, the tar- 
get polarization was changed on a time scale of several 
days, and therefore the corresponding beam polarization 
could be different during these different periods. The 
terms era, a e , a t , and a e t represent the cross sections for 
scattering of polarized electrons off a longitudinally po- 
larized proton in the 15 NH 3 target, while the terms a^ 
and a^ were added to take into account that there are 
additional contributions coming from the 7 unpolarized 
protons bound in the nitrogen nucleus. The asymmetries 
may be written in terms of these quantities as: 
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A et = 
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(15) 
(16) 



where the counts were normalized to the same total 
charge. In these expressions, 



(17) 



and 
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(18) 



The term C-N BG represents the properly normalized 12 C 
and 4 He counts used to simulate and remove the residual 
background contribution 2D 2 {\ + ottj&o m the denom- 
inator. The experimental determination of the normal- 
ization constant C is discussed in the next section. 
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FIG. 5: (Color online) Ml spectra of 15 NH 3 , normalized 12 C, 
and their subtraction for < cos# < 0.2 and invariant masses 
of 1.15 GeV < W < 1.2 GeV (a) and 1.4 GeV < W < 1.5 GeV 
(b). The lines show the 2a cut applied to select -k events. 



region, where one has only events stemming from scatter- 
ing off background nuclei. Since the missing mass spectra 
show a dependence on the center-of-mass polar angle of 
the proton, 6* , the normalization constant was extracted 
for each W and cos 9* bin. Figure \5\ shows a sample of 
missing mass spectra for the 15 NH3, the background, and 
the subtraction. The figure also shows the cut applied to 
select neutral pions. 



E. Multipion background correction at high W 

At higher invariant mass, W, background from the 
lower missing mass tail of multi-pion events contaminates 
the 7T° signal, as can be seen in the right panel of Figurc[5] 
The measured asymmetry, A exp , found by combining the 
number of counts (see Eqs.[l5]and[T6]) is therefore given 
by 



D. Nuclear background subtraction 

The data contains a large background coming from 
scattering off 15 N and the 4 He from the target cooling 
bath. Data on 12 C and 4 He targets were used to re- 
move these contributions. The procedure to determine 
the term CN BG consists of two steps. The first step 
was to determine how to combine the 12 C and 4 He data 
in order to correct for the different amount of 4 He dis- 
placed by the 15 NH3 and 12 C targets. Using a calculation 
based on target thicknesses, densities, and window con- 
tributions, the background counts were constructed as 
n bg = N 12 C _ o.22N 4h °. The second step consists of 
the determination of the normalization constant C. The 
normalization constant is W dependent because of the 
different ratio of protons to neutrons for the 15 N and 12 C 
targets. The constant C was extracted by normalizing 
the missing mass spectra of the 15 N and the constructed 
background to each other in the negative missing mass 
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f^oA^o + f back Aback, 



(19) 



where A^o and Aback are the asymmetries of it and back- 
ground events, respectively, f n o is the ratio of ir° events 
to total events, and fback = 1 — fir° is the fraction of back- 
ground events. Equation [T5] was used to extract A„a . 
Aback w as measured using the same analysis procedure 
used for A exp , except that events were selected with a 
missing mass squared greater than the pion mass squared 
and less than 0.2 GeV 2 . The values of the background 
asymmetry Aback vary from 1% to 30% depending on 
the kinematic bin. The fraction fback was found by fit- 
ting the missing mass spectrum, consisting of the pion 
peak and the background, with two Gaussian curves and 
using those to estimate the number of background and pi- 
ons events in the region selected by the cut. The fraction 
was found to be increasing with the invariant mass as ex- 
pected from a value of 1.2% for 1.25 GeV < W < 1.3 GeV 
up to 12% for 1.6 GeV < W < 1.7 GeV. The overall cor- 
rection fbackAback was from 0.5% to 2.5%. 



P e Pt determination 



where 



The product of the beam and target polarizations, 
P e Pt, was determined by comparing the measured clastic 
asymmetry A meas with the predicted value A t h eo , 
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For the Q 2 range in this analysis, |p- ~ /i p , where fi p 
is the magnetic moment of the proton (/i p = ^g where 
g = 5.5856912). The product of beam and target po- 
larization (P e P t ) was independently determined in six 
Q 2 bins and then the average value of these 6 results 
was calculated. Figure [HI shows the results for the posi- 
tive (P e P?) and negative target polarization data {P e P^). 
These measurements allow one to extract target polariza- 
tions P t a , P 4 b by simply taking the ratio of these products 
and the measured beam polarization P e . 
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FIG. 6: The product |-P e -Pt| as a function of Q for posi- 
tive (filled circles) and negative (open triangles) target po- 
larization runs. The six values for each polarization were 
fitted with a constant in order to obtain the average values 



P e Pt = 0.517 ± 0.002 and P e P t b 



-0.483 ±0.003. 



G. Bin centering corrections 

In order to reduce effects due to the width of the kine- 
matic bins in Q 2 , W, cos 8* , and 4>* , a bin centering cor- 
rection was applied to the data. This correction is model 
dependent and in order to minimize the error, the mod- 



els that describe the data best were chosen. The Sato 
and Lee [J, [35| and MAID07 ,3| models were used for the 
A(1232) and the higher W region, respectively. Neither 
model includes the data presented in this paper in their 
fits. The bin centering correction is defined as: 



/3 



'exact 



(7 



(23) 



where a exact is the model cross section obtained at the 
center values of our bins, and a is the model cross sec- 
tion averaged over the bins in the analysis. This correc- 
tion was calculated for each beam and target polarization 
configuration and applied individually to the number of 
counts iV-ff, -/Vjj, Nfi, and N±± before constructing the 
asymmetries as described in Eqs. 1151 and 1161 



H. Systematic uncertainties 

Several sources of systematic uncertainties were iden- 
tified in the analysis procedure. For most cases, the tech- 
nique used to estimate the size of these uncertainties was 
to vary individual parameters in the analysis, recalculate 
the asymmetries, and compare the new results to the 
original values. 

One of the main sources of uncertainties is due to the 
correction for nitrogen, which was estimated using a 12 C 
target. The statistical uncertainty on the normalization 
procedure is very small, however, the whole procedure is 
based on the assumption that the 12 C target spectrum 
has a similar shape as the 15 NH3 background. Based 
on studies in the elastic region, the normalization factor 
was found to be known with a precision of 3-4%. The 
asymmetry was recalculated using a normalization factor 
changed by 3.5%, and the comparison shows a point by 
point variation of the asymmetry on the order of 4-8%. 

A second source of systematic uncertainties is due to 
the uncertainties of the beam and target polarization val- 
ues. The beam polarization is measured with the M0ller 
polarimcter and it is know with an accuracy of 2%, which 
results in a systematic uncertainty on the asymmetry up 



to 2%. The target polarization, extracted from the prod- 
uct P e P t as described in section IV Fi has a systematic 
and statistical uncertainty associated with it. The sta- 
tistical uncertainty gives a negligible systematic variation 
of the asymmetry, while the systematic error due to the 
uncertainty of the ratio Ge/Gm and to the carbon nor- 
malization gives a total systematic uncertainty on the 
asymmetry of 3.2%. 

Another source of error is the model dependence of the 
bin centering correction. To study this effect a different 
model, MAID03, was used for the correction and the dif- 
ference was assigned as a systematic uncertainty. A point 
by point variation from 1 to 15% was found depending 
on the W range and the asymmetry considered. For the 
most part these uncertainties were negligible compared 
to the statistical error. 

The multi-pion background correction requires an es- 
timate of the number of events contaminating the pion 
sample, fback- The fraction has both a statistical and 
systematic uncertainty that contributes to the systematic 
uncertainty of the asymmetry. The main source of sys- 
tematic uncertainty was the assumption that the back- 
ground does not depend on the polar angle of the pion 
in the center-of-mass, 9* . An overall uncertainty of 20% 
was assigned to fback, to account for fluctuations in dif- 
ferent cos 9* intervals. A change of this amount causes 
a point by point variation in the asymmetry up to 10%, 
most of the time negligible compared to the statistical 
uncertainty. 

The effect of cuts was studied as well, in particular the 
fiducial and missing mass cuts. To study the impact of 
these cuts on the analysis, the asymmetries were mea- 
sured using complementary cuts and compared to each 
other. The two sets of data points were found statisti- 
cally compatible and therefore no systematic uncertainty 
was assigned. 

The overall systematic uncertainties were found by 
combining the different contributions in quadrature and 
are reported point by point in the CLAS Physics Data 
Base |36| . It was found that the systematic uncertainties 
were smaller than the statistical ones for most of the data 
points. 



VI. RESULTS 

Target (At) and double spin (A et ) asymmetries were 
measured as a function of the invariant mass W, the mo- 
mentum transfer squared Q 2 , and the center of mass pion 
angles, 9* and <jf . The data were divided into 8 bins in 
W, 3 bins in Q 2 , 10 bins in cos£?*, and 15 bins in </>*. The 
results consist of 2,435 data points for each asymmetry, 
after eliminating bins with nearly zero acceptance. Ta- 
ble |T] summarizes the binning and Figure [7] illustrates the 
kinematic coverage in Q 2 and W . 

The results were compared with five theoretical ap- 
proaches: MAID07, MAID03, MAID00, and the DMT 
model everywhere, and the Sato and Lee model in the 



Bin 


Range 


Bin size 


# bins 


W 


1.1 - 1.3 GeV 


0.05 GeV 


4 


1.3 - 1.7 GeV 


0.10 GeV 


4 


Q 2 


0.187 - 0.317 GeV" 


0.130 GeV" 


1 


0.317 - 0.452 GeV" 


0.135 GeV" 


1 


0.452 - 0.770 GeV" 


0.318 GeV" 


1 


cos 6* 


-1.0- 1.0 


0.2 


10 


^ 


-180° - 180° 


24° 


15 



TABLE I: Binning for the 1.6 GeV data set. 
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FIG. 7: (Color online) Kinematic coverage in Q 2 and W. The 
dashed lines indicate how the data were subdivided. 



A(1232) region. Figures [8] and [9] show a sample of the 
results as a function of cos 9* . The error bars indicate the 
statistical uncertainties, while the brown boxes indicate 
the systematic uncertainties. The complete numerical re- 
sults are reported in the CLAS Physics Data Base [3a |. 
We see that the models reproduce well the general trends 
of the data. The asymmetry signs and the sign changes 
are generally predicted correctly. Discrepancies are visi- 
ble at the quantitative level and are different for different 
models as discussed in Section IVlAl This shows the sen- 
sitivity of this data set to the specific model ingredients 
for the resonant and non-resonant amplitudes in the first, 
second, and third resonance regions. 

Integrated asymmetries were also extracted as a func- 
tion of the invariant mass W and the momentum trans- 
fer Q 2 . The asymmetries as a function of the invariant 
mass were calculated integrating over the whole range in 
cos#*, Q 2 , and over a partial range in <f>* from 60° to 
156° to avoid regions with acceptance less than 2%. Fig- 
ures [TU] and QT] show the results compared to the models. 
The asymmetries as a function of the momentum trans- 
fer were calculated integrating over the whole range in 
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FIG. 8: The target asymmetries At for Q 2 = 0.385 GeV 2 , 
(j)* = 120°, and invariant masses W = 1.225 GeV (a) and 
W = 1.35 GeV (b). Systematic uncertainties are shown as 
shaded bars. The curves correspond to the MAID07 (solid), 
DMT (dotted), Sato and Lee (dashed), and MAID03 (dashed 
dotted) calculations. 
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tainty. The discrepancy with model parameterizations is 
likely related to an overestimation of the relative strength 
of the dominant A 3 / 2 amplitude in comparison with the 
A 1 / 2 . This sensitivity to the dominant A 3 / 2 amplitudes 
is also shown in Figure 1101 where the magnitude of the 
computed A et increases without the _Di3(1520) resonance 
contribution, bringing it in better agreement with the 
data in the mass range near 1520 MeV. In the mass range 
1.6 < W < 1.7 GeV the Fi 5 (1680) is the dominant res- 
onance. The asymmetry increases rapidly with Q 2 and 
changes sign at Q 2 ~ 0.35 GeV 2 . The MAID07 model ap- 
proximately describes the Q 2 dependence, including the 
sign change. We remark that the observed behavior of 
A et with Q 2 in the region above the A(1232) is consistent 
with the behavior measured in the nir + channel [301 ] . 

The target asymmetry A t depends on the imaginary 
part of interferences involving resonant and non-resonant 
amplitudes, and its interpretation is less straightforward 
than the double spin asymmetry. The asymmetry A t 
is shown in Figure Q2] as a function of Q 2 . There is 
not a simple explanation of the discrepancies between 
the models and data for this observable, as the results 
may depend on the phases of amplitudes that are also 
strongly energy-dependent. In contrast to the double 
spin asymmetry, A t is well described in the mass region of 
the A(1232) only by the SL model, while the DMT and 
MAID parameterizations show significant discrepancies 
with the data. The higher mass regions are reasonably 
well described by the MAID parameterizations, while the 
DMT model shows strong discrepancies, including the 
wrong sign for the four higher mass bins. 



FIG. 9: The target asymmetries A et for Q 2 = 0.252 GeV 2 , 
0*=12O°, and invariant masses W = 1.225 GeV (a) and 
W = 1.65 GeV (b). Systematic uncertainties are shown as 
shaded bars. The curves correspond to the MAID07 (solid), 
DMT (dotted), Sato and Lee (dashed), and MAID03 (dashed 
dotted) calculations. 



cos#*, over a partial range in <f>* from 60° to 156°, and 
over 12 bins in W from 1.1 GeV to 1.7 GeV in steps of 
0.05 GeV. Figures [T2l and [j~3l show the results compared 
to the models. 

Figure [T2l shows the Q 2 -dependence of the double spin 
asymmetry A et for 12 bins in W. The negative asym- 
metry and the approximate linear dependence on Q 2 , 
as well as the negative slope observed in the region of 
the A(1232) resonance are well described by all models. 
In the mass region of the Roper resonance Pn(1440), 
and the L>i3(1520) and S , n(1535) resonances at 1.45 < 
W < 1.60 GeV, the asymmetry is large, positive, and 
increases rapidly with Q 2 . The DMT and MAID models 
all underestimate the magnitude. The strength of the 
Pn(1440) in this low Q 2 range is known to be small p], 
and the 511(1535) has been well measured in 77 electro- 
production. The excitation strength of the Z?i3(1520) 
and the relative strength of the A 1 i 2 and A 3 / 2 ampli- 
tudes of this resonance are the main source of uncer- 
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FIG. 10: Aet as a function of the invariant mass W, integrated 
over the whole range in cos#*, Q 2 , and 60° < (f* < 156°. 
The curves correspond to the MAID07 (solid), DMT (dotted), 
Sato and Lee (dashed), MAID03 (dashed dotted), MAID00 
(long-dashed-dotted) calculations. The long-dashed curve is 
the MAID00 calculation without the Di3(1520) contribution. 
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FIG. 11: At as a function of the invariant mass W, integrated 
over the whole range in cos#*, Q 2 , and 60° < (f>* < 156°. 
The curves correspond to the MAID07 (solid), DMT (dot- 
ted), Sato and Lee (dashed), MAID03 (dashed dotted), and 
MAIDOO (long-dashed-dotted) calculations. 



A. x comparison 

To evaluate in a more rigorous way how well the models 
describe the data, a \ 2 comparison was performed. The 
X 2 was denned as: 
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data 
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mo del \ 2 



data \ 2 



( CT data) 



(24) 



where xf ata is the value of each experimental point for all 
the asymmetries and xf 10 ^ 1 is the corresponding value of 
the theoretical prediction. Since the model is given with- 
out errors, only the experimental statistical uncertainties 
<7 data were used in the denominator. The comparison 
gave the results listed in Table [TTJ The calculation led us 



Model 


At 


A et 


A t 


Aet 




W < 1.3 GcV 


W > 1.3 GcV 




ndf = 1306 


ndf = 1129 


MAID07 


1.98 


1.04 


0.99 


1.30 


MAID03 


1.75 


1.05 


1.20 


1.86 


MAIDOO 


1.97 


1.03 


1.97 


1.07 


SL 


1.00 


1.05 


- 


- 


DMT 


2.09 


1.02 


1.71 


1.17 



TABLE II: x P er number of degree of freedom (ndf) com- 
parison between the data and the five models. 

to make a few conclusions about the models. First of all, 
in the A(1232) region, all the models are in good agree- 
ment with the measured double spin asymmetry, A et . 



This is not surprising since A et is dominated by the res- 
onance multipoles (i.e. \M 1+ \ 2 , Re{EJ + M 1+ }) [H[ that 
dominate the unpolarized cross section. Since the model 
fits are based on unpolarized cross sections, this agree- 
ment is expected. At higher W none of the models give 
a consistently good description of the data for both the 
target and the double spin asymmetry. The MAID07 pa- 
rameterization shows improvements for A t over previous 
versions, while A et is not as well reproduced as by the 
older MAIDOO version. At this level it is not possible to 
draw further conclusions and only including these data 
in a new global fit for these models will give a better 
understanding of the sensitivity of the data and possibly 
identify the origin of the discrepancies. 



VII. SUMMARY 

In this paper we have presented a set of longitudi- 
nal target polarization asymmetries A t and beam-target 
double polarization asymmetries A et for the reaction 
ep — > epir covering the first, second, and part of the 
third nucleon resonance regions in a wide range of az- 
imuthal and polar center-of-mass angles for the final state 
7T°. The data also cover a range in photon virtuality 
Q 2 . For the first time, target polarization asymmetries 
have been measured in the resonance region and in the 
mass range above the A(1232) resonance. The new data 
complement the alr eady published results obtained in 
the A(1232) region [13| for different kinematics. Lon- 
gitudinal target asymmetries and beam-target asymme- 
tries were measured in 2,435 bins each. The polarization 
asymmetries show strong sensitivity to both resonant 
and non-resonant amplitudes and are highly selective of 
model paramctcrizations for amplitudes underlying the 
7T° electroproduction process. Due to the large number 
of data points, only samples of specific kinematics and 
integrated quantities are compared with model predic- 
tions. Although the models agree well with the double 
spin asymmetries in the A(1232) mass region, three of 
the four models surveyed have a substantial disagreement 
beyond the A(1232) region and all models show a signifi- 
cant disagreement with the single spin target asymmetry 
data. The discrepancies, when quantified using the sta- 
tistical interpretation from Ref. [37| of the x 2 values in 
Table |TT] above, amount to less than a 5% probability 
that the measured double spin asymmetry data set (A et ) 
would be consistent with any of the given models for 
W above the A(1232). For the single spin target asym- 
metry (A t ) the probability is less than 59%. The full 
impact of the complete data set on our understanding of 
the nucleon structure in the regime of strong QCD and 
confinement can only be obtained by including these new 
data sets in global analyses that incorporate all exclusive 
pion electroproduction cross sections and polarization ob- 
scrvables. This effort is currently underway, and is the 
subject of forthcoming publications [38[. The full set of 
data is available at the CLAS Physics Data Base [361 ] . 
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FIG. 12: A e t as a function of the photon virtuality Q 2 , integrated over the whole range in cos #*, 60° < (j>* < 156°, and 12 
invariant mass regions from 1.1 GeV to 1.7 GeV in steps of 0.05 GeV. The curves correspond to the MAID07 (solid), DMT 
(dotted), Sato and Lee (dashed), and MAID03 (dashed dotted) calculations. 
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FIG. 13: A t as a function of the photon virtuality Q 2 , integrated over the whole range in cos#*, 60° < (f>* < 156°, and 12 
invariant mass regions from 1.1 GeV to 1.7 GeV in steps of 0.05 GeV. The curves correspond to the MAID07 (solid), DMT 
(dotted), Sato and Lee (dashed), and MAID03 (dashed dotted) calculations. 
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